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ABSTRACT 

We have measured star formation histories (SFHs) and stellar masses of galaxies 
detected by the Balloon-borne Large Aperture Sub-millimetre Telescope (BLAST) 
over ~ 9 deg 2 centred on the Chandra Deep Field South. We have applied the recently 
developed SFH reconstruction method of Dye et al. to optical, near-infrared and mid- 
infrared photometry of 92 BLAST galaxies. We find significant differences between the 
SFHs of low mass (<J10 1:L M ) and high mass (^lO 11 !^) systems. On average, low 
mass systems exhibit a dominant late burst of star formation which creates a large 
fraction of their stellar mass. Conversely, high mass systems tend to have a significant 
amount of stellar mass that formed much earlier. We also find that the high mass SFHs 
evolve more strongly than the low mass SFHs. These findings are consistent with the 
phenomenon of downsizing observed in optically selected samples of galaxies. 
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1 INTRODUCTION 

Encoded in every galaxy's spectrum is a record of its en- 
tire life from birth, up to the epoch at which it is observed. 
Unlocking this information, by determining the variation in 
star formation rate (SFR) with age, is a key step towards 
understanding how galaxies form and evolve. The measure- 
ment of star formation histories (SFHs) therefore plays a 
crucial role in the development of an accurate model to de- 
scribe the range of processes experienced by galaxies and 
the subsequent formation of stellar mass in the Universe. 

Currently, the most elusive population of galaxies are 
systems heavily obscured by dust, often detected only at 
sub-millimetre (submm) wavelengths. Approximately half 
of all light emitted by galaxies is absorbed by dust and 
re-radiated in the submm. However, compared with studies 
at optical wavelengths, little is known about submm galax- 
ies, in particular, how they relate to local systems. Submm 
selected samples of galaxies therefore provide an unavoid- 
ably important set of constraints on a complete and self- 
consi stent view of gal axy formation mechanisms. 

(|Dve et al. I l2008h conducted an analysis of the SFHs 
of 850 fim selected galaxies and found that these systems 
are typically dominated by a strong burst of star formation 
late in their history, but that around half of their stellar 
population had already formed over the first half of their 
lives. This study also revealed a surprising deficit of high 
mass (> 5 x lO n M0) systems at redshifts z < 2, strong 
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evidence of downsizing in the population, possibly explained 
by the evolution of these systems into massive ellipticals. 

The Balloon-borne Large Aperture Sub-millimetre Tele- 
scope (BLAST) recently provided a catalogue of hundreds of 
subm m selected galaxies over ~ 9 deg 2 of sky l|Devlin et al. I 
l200Sh centred on the Chandra Deep Field South (CDFS). 
The identification of radio and 24 /j ,m counterparts to the 
majority of sources (iDve et al. Il2009 | ) combi ned with follow- 
up optical spectroscopy ( Eales et al. 1120091 ) and a compre- 
hensive suite of archival optical, near-infrared (near-IR) and 
mid-infrared (mid-IR) imaging results in these data being 
the largest, most thoroughly characterised and carefully pro- 
cessed sample of 250-500 /J.m selected sources to date. The 
sample and its supporting multi-wavelength data therefore 
presents a perfect opportunity to conduct a study of SFHs 
of submm selected galaxies. 

The purpose of this letter is therefore to carry out an in- 
vestigation of the SFHs and stellar masses of BLAST galax- 
ies in a simil ar vein to the stud y of 850 fim selected galaxies 
conducted bv lDve et al. I l|200ct ). To c ompute SFH s, we have 
used the recently developed method o f |Pvell|200Sl ). We have 
applied the method to optical, near-IR and mid-IR photom- 
etry of the counte rparts to the BLAST sources identified by 
IDve et al. I \200$ ). 

In Section [2] we briefly outline the procedure used for 
applying the SFH reconstruction method. Section [3] de- 
scribes the data. The results are presented in Section [4] 
Finally, we summarise in Section [5] Throughout this let- 
ter, the following cosmological parameters are assumed; 
H = 71 kms -1 Mpc -1 , ft m = 0.27, ft A = 0.73. 
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2 SFH RECONSTRUCTION METHOD 

I Dye I l|2008l ) provides a detailed description of the SFH re- 
construction method. The purpose of the brief outline pre- 
sented here is both for completeness and to give specific 
details of the procedure we have used in this implementa- 
tion. 

2.1 Method Outline 

The method divides a galaxy's history into discrete blocks 
of time. This results in a relatively low resolution SFH, but 
one that does not adhere to a prescribed (i.e., potentially 
biased) parametric form. Using a synthetic library of simple 
stellar population (SSP) SEDs, the fluxes resulting from a 
constant SFR normalised to one solar mass in each block, as 
measured in the observer frame across a range of filters, are 
calculated. Finding the contribution of flux from each block 
in each filter that best fits a set of observed fluxes is a linear 
problem that can be solved exactly. The solution directly 
yields the galaxy's SFH and stellar mass. In th i s lette r, we 
have used the stellar SED libraries of lMarastonI (|2005r i and, 
for comparison, iBruzual fc Chariot I (|2003h . 

Starting with a SSP SED, if p , of metallicity Z, a com- 
posite stellar population (CSP) SED, LI, is generated for 
the ith block of constant star formation in a given galaxy 
using 



LI = 



1 

Ah 



dt'Ll sp (r(z)-t') 



(1) 



where the block spans the period to U in the galaxy's 
history and r is the age of the galaxy (i.e., the age of the 
Universe today minus the look-back time to the galaxy). 
The quantity Ati = ti — t%-i ensures that the CSP is nor- 
malised to one solar mass. The above integral is evaluated 
by interpolating linearly in log(t) between the discrete time 
intervals at which the SSP SEDs are given in the libraries. 
Note that the method assumes that Z does not vary with 
age. 

To model the effects of extinction on the final SED (i.e., 
the SED from all blocks in the SFH), reddening is applied. 
This is achieved by individually reddening the CSP of each 
block using L\ >R = h{ i(r - 4fc ( A > A v/«v w here, A v is the 
extin ction, Rv = 4.05 and the Calzetti Law l|Calzetti et al. I 
2000) is used for fc(A). The model flux (i.e., photon count) 
observed in filter j from a given block i in the SFH when 
the galaxy lies at a redshift z is then 
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where is the luminosity distance and Tj is the transmis- 
sion curve of filter j. 

To find the normalisations which result in a set of 
model fluxes that best fits the observed fluxes, the following 
X 2 function is minimised 



x 2 = E- 



O'i-F'ij 



F° 



(3) 



where F° bs is the galaxy flux observed in filter j and cjj 
is its error. The sum in i acts over all A'biock SFH blocks. 
The minimum x 2 occurs when the condition dx 2 /da,i — 



is simultaneously satisfied for all a^. This gives a set of 
equations linear in the dj which ar e solvable using a standard 
matrix inversion (see I Dye 1 120081 . for more details) . The a» 
are the stellar masses formed in each block so that the total 
stellar mass of the galaxy is the sum M* = ^ JVblock ai _ The 
SFRs and hence the SFH is then given directly by dividing 
the di by the time spanned by each corresponding block. 
Formal errors on the a,i are obtained from the covariance 
matrix, computed in a simple additional step. To allow for 
uncertainty in source redshift, we performed a Monte Carlo 
simulation, randomising the redshift according to its error, 
and combined the resulting scatter in the ai in quadrature 
with the formal errors. The total error on the stellar mass 
was obtained in the same mann er. 

As discussed in I Dye I l|2008l ). regularisation must be ap- 
plied to ensure that the linear solution is well defined. The 
strength of regularisation is controlled by a parameter re- 
ferred to as the regularisation weight, denoted w hereafter. 



2.2 Fitting Procedure 

The procedure outlined in the previous section is a single 
linear step which computes the SFH that gives the best fit 
(i.e., minimum x 2 ) to an observed set of fluxes for a given set 
of parameters z, Z, Av, A'biock and w. This step is nested 
inside a non-linear search for the set of parameter values 
that gives the best ov erall fit to the observed fluxes. As dis- 
cussed in IDve I ()2008l ). finding the best global solution can 
not be achieved by minimising x 2 because the effective num- 
ber of degrees of freedom depends on w in an unquantifiable 
manner. It is therefore not possible to use x 2 to make a fair 
comparison of the goodness-of-fit between two parameter 
sets with differing values of w. To make a fair comparison, 
one must turn to Bayesian statistics and treat regularisa- 
tion as a prior. In this way, the Bayesian evidence, denoted 
e hereafter, allows different sets of parameters to be ranked 
fairly. The best global solution is that which maximises e. 

We investigated a range of different schemes to max- 
imise e. The most efficient and reliable scheme that we found 
combines a standard grid search with a downhill simplex 
minimisation to find the minimum of — In e (hence the max- 
imum of e). We use the simplex routine, linearly computing 
the SFH and evaluating e each time, to find the best pair of 
values of Av and logui whilst keeping Z and A'biock fixed. 
Then, at the outer-most level, we step through a grid of reg- 
ularly spaced trial values of log Z and A'biock over the ranges 
-2 < \og 10 {Z/Z Q ) < 0.3 and 2 < A'biock < 5. The value of z 
was fixed at the measured redshift of the galaxy (see Section 
at all times. 



3 DATA 



I Dye et ai~l (|2009T ) found optical counterparts to 114 of the 
~ 130 radio- and/or 24 /xm-identified BLAST sources lo- 
cated in the region covered by their optical data. All 114 
sources were detected at 5<r in at least one of the three 
BLAST bands (250, 350 and 500 fim). As described below, 
we acquired a range of optical, near-IR and mid-IR pho- 
tometry for these 114 sources. We rejected 12 sources on 
the basis of having less than five photometric data points, 
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applying our analysis to the remaining 102. In all cases, we 
used total magnitudes/fluxes. 

Optical photometry was taken from either the 
Spitzer Wide-area In fraRed Extragalactic survey (SWIRE; 
lLonsdale et~ai~ll2004l ), or the 17 band COMBO-17 (Classify- 
ing O bjects by Mediu m-Band Observations in 17 filters) sur- 
vey l|Wolf et al. 1120041 ). The SWIRE survey covers ~ 5deg 2 
to a depth of r ~ 24.5 (Vega, 5<r) and COMBO-17 covers 
~ 0.25 deg 2 to R ~ 26.0 (Vega, 5<r). The SWIRE optical 
catalogue directly provides total magnitudes. For COMBO- 
17, we computed total magnitudes from aperture fluxes in 
all filters using the given source-specific correction derived 
in the R band. For objects detected in both surveys, we 
amalgamated both sets of photometry. 

For the near-IR, we used the J and K photometry 
from the Multi-wave length Survey by Yale-Chile (MUSYC 
iGawiser et ai~ll2006l) . This reaches depths of 22.1 and 20.5 
(Vega, 5cr) in J and K respectively. Although MUSYC 
only covers ~ 0.3 deg 2 , the survey area is centred on the 
CDFS where the deep BLAST observations were made (see 
iDve et al. 120091) . hence we obtained near-IR photometry for 
52 of the 102 sources. The MUSYC catalogue gives aperture 
and total K band fluxes. We obtained total J band fluxes 
by scaling the aperture fluxes by the totakaperture K band 
flux ratio. 

Mid-IR photometry was taken from the SWIRE sur- 
vey which completely encompasses the region containing the 
BLAST sources considered in the present work. We used 3.6 
and 4.5 /^m fluxes which are limited to a 10a sensitivity of 
10 /iJy at 3.6 /im and a 5<r sensitivity of 10 /^Jy at 4.5 /im. 
For reconstructing source SFHs, we only used 3.6 fim pho- 
tometry for objects with redshifts z > 0.6 and additionally 
4.5 fim photometry only for objects at z > 1. These redshift 
limits ensure that the observed photometry does not extend 
beyond the restframe K band where the empirical stellar 
component of the Maraston SEDs ends. 

We imposed a minimum photometric error of 0.05 mag 
for all photometry to allow for zero point uncertainties and 
mismatches in calibration between different datasets. 

In terms of redshifts, 47 of the sources have spec- 
troscopic redshifts taken from iMoncelsi et al. I l|201Ch . For 
the rema inder, we used the photometric redshifts as- 
signed in IDve et al. I dgOOSh. These are taken from either 



iRowan-Robinson et al. I (|2008l ) which uses the SWIRE opti- 
cal photometry as well as the SWIRE 3.6 and 4.5 /im data 
or from COMBO-17. 



4 RESULTS 

We applied the procedure outlined in Section 12.21 to max- 
imise the Bayesian evidence for each galaxy. We found that 
5-10% of reconstructions resulted in very poor fits to the 
observed photometry and/or gave unphysical SFHs (i.e., 
strongly negative bursts of star formation). These failures 
were readily discarded by applying a cut in the evidence of 
In e > —100. This removed 10 of the 102 sources, includ- 
ing three low redshift galaxies (z < 0.1) which are optically 
well resolved (up to ~ 50" in diameter) and therefore sus- 
ceptible to strong photometric biases between the different 
catalogues. All analysis that follows in this letter is applied 
to the remaining 92 systems. 
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Figure 1. Example rest-frame SEDs fitted to optical and near- 
IR photometry of two BLAST sources. The darker SED, corre- 
sponding to a z ~ 0.8 galaxy with observed photometry plotted 
as filled photometric data points, is fit with an SFH where 60% 
of its stellar mass was formed in the last 10% of its history and 
15% formed in the first half. The lighter SED, corresponding to 
a 2 ~ 1.2 galaxy with observed photometry plotted as unfilled 
data points, has a SFH with only 3% of its stellar mass formed 
in the last 10% of its history and 85% formed in the first half. 
Fluxes are in arbitrary units and SEDs are independently scaled 
for clarity. Inset panels show the SFR versus fractional galaxy age 
from birth (t = 0) to the epoch its observed (t = 1) for each SED. 



Our findings indicate that the vast majority of BLAST 
galaxies have SFHs which peak at late times. Around 10% 
of systems are more consistent with early type SEDs having 
formed at least three quarters of their stellar mass in the 
first half of their lives and less than 2% of their mass in the 
last tenth of their lives. Figure [1] shows two fairly extreme 
example SEDs of BLAST sources, one dominated by late 
star formation and the other having formed nearly all of its 
stars early on, to illustrate the range SEDs observed. 

Figure [5] shows how the stellar masses vary with red- 
shift. There is an obvious trend of increasing mass with in- 
creasing redshift, limited by the rarity of high mass galaxies 
at one end of the scale and the sensitivity to low luminosity 
at the other. The few low mass outliers are systems domi- 
nated by extremely late and intense bursts of star formation 
and hence have low mass-to-light ratios. Four of the 92 ob- 
jects are flagg ed as having a dominan t active galactic nuclear 
component bv lMoncelsi et al. I l|2010l ) based on their spectral 
line strength ratio [NII]658.3/Hq. These four objects are all 
located near the high mass edge of the envelope shown in 
Figure [2] and spread throughout the redshift range. 

We investigated the possibility of evolution in the SFHs 
and whether this depends on mass by segregating the sample 
of 92 galaxies into four approximately equally sized sub- 
samples divided at z = 0.5 and M, = 9 x 10 10 M Q . Figure 
[3] shows the average SFH for each sub-sample, rebinned to 
a common resolution of five blocks. In the averaging, we 
normalised the SFH of each galaxy to units of fractional total 
stellar mass formed per fractional age so that the integral of 
SFH over fractional age (i.e., the total stellar mass) is equal 
to unity. 
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Figure 2. Stellar mass in units 10 Mg versus redshift for the 
BLAST sources. The dashed lines at z = 0.5 and 9 X 10 10 M di- 
vide the sample into four approximately equal sub-samples which 
we have then used in our analysis of SFH evolution by mass. 
Open points indicate sources where a spectroscopic redshift has 
been used and filled points with errors show the remaining sources 
with photometric redshifts. Mass errors include photometric and 
redshift uncertainties. 
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Figure 3. Mean SFHs of the four sub-samples of BLAST galaxies 
as delineated by the redshift and stellar mass limits indicated (see 
Figure [2{. The error bars on each histogram bin indicate the la 
error on the mean which includes the formal error from the linear 
inversion and the uncertainty in source redshift computed in the 
Monte Carlo analysis (see Section |2. 111 . All SFRs are expressed in 
units of fraction of total stellar mass per fractional age (i.e., the 
area under each histogram in these units is 1). The fractional age 
varies from t = 0atz = ootot = lat the epoch of the galaxy's 
redshift. 



The most striking feature seen in the plots is the differ- 
ence between the low and high mass sub-samples. On aver- 
age, both low mass SFHs are dominated by a late burst of 
strong star formation activity accounting for ~ 40% of the 
total stellar mass. Conversely, both high mass SFHs show 
that a much smaller fraction of stellar mass (~ 5 — 10%) is 
created during this last period, the majority of mass being 
formed at earlier times. 

Another obvious effect seen in Figure [3] is that the high 
mass sources exhibit a more prominent difference in their 
SFHs in moving from high to low redshifts than the low 
mass sources. The high mass sources have therefore, by this 
definition, undergone more evolution. To quantify the signif- 
icance of this, we computed the reduced \ 2 statistic between 
the low and high redshift SFHs for the low and high mass 
sub-samples in turn. For the high mass SFHs, the statistic 
is Xr — 3.57 ±0.63 compared to Xr = 0.51 ±0.63 for the low 
mass SFHs. The change at high mass is therefore significant 
at the ~ 3<J level whereas the low mass source SFHs are 
consistent with no change. This is synonymous with down- 
sizing where the instantaneous star formation rate in high 
mass galaxies evolves more strong ly than that in low mass 
systems (e.g.. iHeavens et al. 1120041 ). 

The SFHs computed in terms of fractional mass and 
fractional galaxy age are a very useful diagnostic since they 
effectively normalise out the large scatter in mass and red- 
shift present in the necessarily coarsely binned sub-samples. 
This makes the mean trends more conspicuous. However, to 
compare with more traditional studies of the evolution of 
star formation, we estimated instantaneous absolute SFRs. 
For each source, we computed a 'pseudo-instantaneous' SFR 
by dividing the absolute stellar mass created in the last 
SFH block by the real time spanned by the block. We 
found that the pseudo-instantaneous SFR for the high mass 
sources changed from 75 ± 26 Mq yr _ at high redshifts to 
20 ± 5M©yr~ 1 at low redshifts. In comparison, the change 



for the low mass sources is from from 43 ± 23 Mq yr _1 at 
high redshifts to 9 ± 2 Mq yr~ at low redshifts. The con- 
clusion is therefore that we detect no significant difference 
in the evolution of the pseudo-instantaneous SFR between 
the high and low mass sources. To detect an absolute trend 
such as this, more sources would be required to enable finer 
binning in mass and redshift. 

An interesting point to note is that the rate of forma- 
tion of stellar mass, which is highest at early and at late 
times in the high mass, hig h redshift sub-sam ple, is very 
similar to that measured by iDve et al. I HH) for 850 /im 
selected sources. This is perhaps not too surprising given 
the large overlap of this sub-sample with the 850 fim sample 
which has a median value of redshift and log 10 (M/MQ) of 
1.6 ± TO and 11.5 ± 0.5 respectively, where the errors give 
the standard deviation. The fact that such a large stellar 
population was already in place at higher redshifts suggests 
that the peak star formation rate occurred significantly ear- 
lier in the history of the Universe for high mass systems 
than for low mass sys tems. This behaviour was observed by 
I Heavens et al. I (|2004l ) for optically selected galaxies. 

To verify the robustness of our results, we conducted 
a series of tests. The first was to see if the inferred SFHs 
are intrinsic or merely the effect of reddening. For exam- 
ple, an intrinsically late-type galaxy with strong reddening 
could give rise to a reconstructed SFH with artificially sup- 
pressed late star formation. We therefore plotted the frac- 
tion of mass formed in the last 10% of each galaxy's history, 
Mio%, against Ay- Since late activity strongly dominates 
the shape of the observed SED, M 10 % is a sensitive indica- 
tor of SED type. Therefore, a strong degeneracy between the 
inferred lateness of an SED and extinction would manifest 
itself as an obvious positive correlation between M 10 % and 
Av (if the SED is more reddened, more late stellar mass 
is required to maintain a fit to the observed photometry). 
Figure [4] shows these two quantities plotted for all galaxies 
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Figure 4. The relationship between the fraction of mass formed 
in the last 10% of each galaxy's history, M 10 %, and Ay- 



in our sample. The scatter is large, although there is some 
evidence of a correlation. For example, all galaxies which 
form less than 10% of mass in the last 10% of their history 
have values of Ay < 0.5. However, this is at least partly 
explained by the simple fact that elliptical galaxies tend to 
have little or no dust. We therefore conclude that any degen- 
eracies between reddening and reconstructed SFHs do not 
significantly affect the results described in this letter. 

The second test addresses the concern that there may 
be a potential bias introduced by including the mid-IR pho- 
tometry based on redshift. We therefore isolated all sources 
whose SED fitting used mid-IR photometry and repeated 
the fitting without it. Within the errors, which were made 
larger by the lack of mid-IR photometry, we found negligible 
differences in the results. 

As a third t est, we repeated the analys is using the stellar 
SED library of ijBruzual fc Chariot lliooj ). The only signifi- 
cant difference we found was that the stellar masses were an 
average of 40% higher than those computed with the Maras- 
ton SEDs. The reconstructed SFHs of the high mass sources 
also showed slightly higher fractions of mass at early times, 
but this was within the errors. 

Our final test was to repeat the analysis without the 
four sources flagged as having strong AGN activity. The dif- 
ferences were again negligible. 



5 SUMMARY 

Using optical, near-IR and mid-IR photometry, we have re- 
constructed the SFHs of a sample of 92 submm sources de- 
tected by BLAST. Their SFRs peak at late times on average, 
consistent with the high instantaneous SFR inferred from 
their submm emission. We divided the sample by mass and 
redshift into four sub-samples. Our findings clearly indicate 
that the low mass sources form a much higher fraction of 
their stellar mass at late times than the high mass sources, 
consistent with the notion that the SFRs of higher mass 
galaxies peaked at earlier times. Furthermore, the SFHs of 
the higher mass sources evolve more strongly than the SFHs 
of lower mass sources. This behaviour is synonymous with 
downsizing observed in optically-selected samples of galax- 
ies. Finally, our high mass, high redshift sub-sample shows 
evidence of stellar mass being formed predominantly at late 
and at early times, but less so when the galaxies are middle- 



aged. The same trend wa s also observed in the sample of 
850 /mi selected sources bv lDve et al. 1 ((2008), although this 
is perhaps not surprising given the similar range of masses 
and redshifts in the sub-sample. 

This letter has analysed the optical counterparts to ap- 
proximately one third of the full BLAST detected sample 
of galaxies presented in iDve et al~1 ([2009) . It is now un- 
derstood that a significant fraction of the galaxies detected 
solely at 500 /im, are t he result of flux boosti ng and therefore 
probably not real (see lMoncelsi et al. | [2010). Our sample of 
92 galaxies contains five 500/im-only detections. Based on 
the findings of Moncelsi et al., we expect that only ~ 1 of 
these is not real. Of the BLAST sources believed to be real, 
around 50% of these were not identified with radio and/or 
24 /im counterparts. Whilst some of these may have been 
detected in the optical/near-IR surveys used in this letter, 
it is likely that the majority will be heavily dust obs cured 
systems lying at higher redshifts fsee lDve et al. |[2009l ). 

Although we have allowed for extinction by dust, 
there could be regions in the BLAST galaxies com- 
pletely obs c ured a t rest-frame optical/near-IR wavelengths. 
IDve et al~l (|2008l ) made a correction for this effect using 
far-IR/submm bolometric luminosity, finding that fully ob- 
scured star formation could result in the generation of up 
to an additional 50% of the galaxy's total stellar mass. It 
is likely that this fraction is even larger for those BLAST 
galaxies with radio/24 /im counterparts but no optical/IR 
counterparts. 

The findings described in this letter represent a taster 
of what would be possible with a significantly larger sam- 
ple of sources. Although the BLAST data currently offer 
the largest and most thoroughly processed sample of galax- 
ies selected over the wavelength range of 200 to 600 /mi to 
date, the increased sensitivity and resolution of the Her- 
schel space observatory, which recently started operation, 
will soon provide vastly increased numbers of sources. This 
will enable significantly reduced uncertainties and therefore 
much improved constraints on models of galaxy evolution 
and formation. Nevertheless, the BLAST data will still pro- 
vide a very valuable benchmark for the Herschel data and 
the various analyses that will emerge for some time to come. 
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